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J. Graf von der Palilen and D. Tsiklauri, Phys. Plas. 21, 060705 (2014), established that the 
generation of octupolar out-of-plane magnetic field structure in a stressed X-point collapse is due 
to ion currents. The field has a central region, comprising of the well-known qaudrupolar field 
(quadrupolar components), as well as four additional poles of reversed polarity closer to the corners 
of the domain (octupolar components). In this extended work, the dependence of the octupolar 
structure on domain size and ion mass variation is investigated. Simulations show that the strength 
and spatial structure of the generated octupolar magnetic field is independent of ion to electron 
mass ratio. Thus showing that ion currents play a significant role in out-of-plane magnetic structure 
generation in physically realistic scenarios. Simulations of different system sizes show that the 
width of the octupolar structure remains the same and has a spacial extent of the order of the ion 
inertial length. The width of the structure thus appears to be independent on boundary condition 
effects. The length of the octupolar structure however increases for greater domain sizes, prescribed 
by the external system size. This was found to be a consequence of the structure of the in-plane 
magnetic field in the outflow region halting the particle flow and thus terminating the in-plane 
currents that generate the out-of-plane field. The generation of octupolar magnetic field structure 
is also established in a tearing-mode reconnection scenario. The differences in the generation of the 
octupolar field and resulting qualitative differences between X-point collapse and tearing-mode are 
discussed. 

PACS numbers: 52.65.Rr;52.30.Cv;52.27.Ny;52.35.Vd;52.35.Py 


I. INTRODUCTION 

Magnetic Hall reconnection, first proposed by B. Son- 
nerup [l| , is a mode of reconnection relying on the decou¬ 
pling of ions and electrons in a diffusion region and is of 
great interest in the study of magnetic reconnection. It 
presents an alternative to the the Petschek model, which 
relies on an anomalous resistivity [2j. Even in setups 
suitable for Petschek reconnection, contributions of Hall 
effects need to be considered. A recent analytical result, 
corroborated by a numerical study, shows that the transi¬ 
tion from Petschek to Hall reconnection occurs when the 
half-length of the current sheet reaches the ion inertial 
length [3| . This was shown to be a direct consequence of 
a generalised scaling law, relating the reconnection rate 
to the distance between the X-point and the start of slow 
mode shocks. 

An observational consequence of Hall reconnection is 
the generated quadrupolar out-of-plane magnetic field, 
induced by currents resulting from the decoupling of elec¬ 
trons from ions i.e. the Hall currents, first demonstrated 
in a study by Teresawa [ 4 ] . The effect was further shown 
to occur in numerical Hybrid simulations and later 
in a full Particle In Cell (PIC) numerical simulation [gj. 
However, as shown in Ref. [§], kinetic simulations of mag¬ 
netic reconnection, where the Hall term was excluded, 
can also lead to quadrupolar magnetic field structure 
generation, due to ion diamagnetic drifts driven by an 
anisotropic ion stress tensor. By being an observational 
signature of magnetic reconnection, the quadrupolar field 
has thus been of great interest in recent spacecraft mis¬ 
sions, including Polar 0 and Cluster 0]. Both mis¬ 


sions observed individual magnetic poles in the magneto¬ 
tail of the Earth. Subsequently, a full quadrupolar pat¬ 
tern was observed in a multi-spacecraft Cluster mission 
0 . The experimental evidence of the full quadrupolar 
structure was also found at the MRX facility [ijj]. How¬ 
ever, it was shown in Ref. 0 that inhomogeneous ion 
flow and pre-existing out-of-plane magnetic fields, can 
lead to the generation of quadrupolar out-of-plane mag¬ 
netic field structure without the Hall term. Thus, the 
generation of a quadrupolar magnetic field is not neces¬ 
sarily a tell-tale sign of Hall-mediated magnetic recon¬ 
nection. 

Ref. fl5l ] proposes an analytical model, explaining the 
Hall out-of-plane quadrupolar magnetic field near an X- 
point as the result of electron motion towards and away 
from the X-point, as field lines reconnect. A uniform ion 
distribution was assumed. Since ions decouple from the 
magnetic field sooner than electrons, they move indepen¬ 
dently of field lines near the X-point. On the other hand, 
electrons are assumed to be coupled to the field lines and 
thus only move with the field and along the field lines. 
Since the spacing of field lines increases as they approach 
the X-point, the electron density decreases. Thus, due 
to the uniform ion density, this results in a net positive 
charge. Electrons in the inflow region therefore move 
along the field lines to towards the X-point to restore 
charge neutrality and then move away from the X-point 
in the out flow region, leading to a quadrupolar pattern. 

While their model describes the generation of the 
quadrupolar field, Ref. [l5j points out that this model is 
limited by the non-inclusion of ion currents, resulting in 
a quadrupolar field that stretches along separatrix arms 
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FIG. 1: (Left panel) Reconnection at an X-point super¬ 
imposed with electron motion (as indicated by the labelled 
track) and the resulting out-of-plane magnetic field structure, 
as given by the analytical model in Ref. 0- Black arrows 
on field lines signify inflow and outflow regions. Ions here 
are assumed to be decoupled and uniformly distributed. As 
the spacing of field lines increases at the X-point, coupled 
electrons in the inflow region move towards the X-point to 
restore charge balance. Similarly, the electrons in the outflow 
region move outwards along the field lines and the characteris¬ 
tic quadrupolar out-of-plane magnetic field structure is gener¬ 
ated. Due to the non-inclusion of ions, out-of-plane magnetic 
structure is not localised and extends along the seperatrices. 
(Right panel) The setup and resulting current and out-of¬ 
plane magnetic field generation for X'-point collapse, from 
simulations in Ref. 16|. It is shown that, as in-plane field 
lines reconnect, electrons move towards and then away from 
the X-point, generating quadrupolar structure as described in 
Ref. 15]. Away from the X-point, ions move independently 
from the field (see dashed tracks), generating a magnetic field 
of opposite polarity to that of the quadrupole, thus gener¬ 
ating an overall octupolar structure. Electrons at the edge 
of the ion diffusion region move such that they cancel the 
out-of-plane field, thus making it localised. 


indefinitely. In a recent simulation study by J. Graf von 
der Pahlen and D. Tsiklauri 0 it is shown that, in an 
X-point collapse scenario, ion currents not only provide 
a cut-off to the quadrupolar field, but contribute to the 
out-of-plane magnetic structure themselves. The out-of- 
plane magnetic field that emerged in Ref. [Tfit was shown 
to have the well-known quadrupolar field at the centre, 
generated by electron currents, and four regions of op¬ 
posite magnetic polarity on the outside, resulting from 
ion currents, as illustrated in Fig. [I] The overall field 
appears to have an octupolar structure (we shall refer 
to the inner magnetic quadrupolar field as quadrupolar 
components and the outer field of opposite polarity as 
the octupolar components). 

At the beginning of the simulation, a different type 
of octupolar magnetic field, with smaller field strength 
(~ 3% of external in-plane field), emerges (see Fig. 2 
in Ref. 0 ])- The same effect was demonstrated in 
Ref. l!| (see their Fig. 5) using a hybrid simulation. 
However, the present study focusses on the larger field 
structure (« 15% of external in-plane field), emerging 
later in the simulation. Further, a octupolar signatures 
have been observed in tearing-mode reconnection scenar¬ 


ios with multiple islands as shown in Ref. 0, Fig. 4, 
and in Ref. [13]), Fig. 3. The emergence octupolar struc¬ 
ture in these scenarios is linked to the island coalescence 
and is beyond the scope of this paper. 

Here, we extend Ref. 0 ] by investigating the depen¬ 
dence of the octupolar field on variation of electron to ion 
mass ratio and the domain size. It is established that the 
generation of octupolar field is neither the result of un¬ 
physical boundary conditions nor unrealistic mass ratios. 
Thus we show that the results are relevant for real labo¬ 
ratory experiments and spacecraft observations. Further, 
by simulating a tearing-mode set-up, as previously stud¬ 
ied in Ref. !8j, it is shown that a similar type of octupo¬ 
lar structure can also be found in reconnection scenarios 
other than A'-point collapse, where their emergence have 
been previously overlooked. 


II. SIMULATION SETUP 

Previous works on collisionless X-point collapse can be 
found in Refs. 0,0-0. In the simulation results pre¬ 
sented in this study, the in-plane magnetic field is that 
of a standard X-point collapse configuration, first intro¬ 
duced by Dungey in 1953 [22|, given by 
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where Bq is characteristic magnetic field strength, L is 
the global external length-scale of reconnection, and a is 
the stress parameter (see e.g. chapter 2.1 in Ref. 0 )- 
A uniform current is imposed at time t = 0 in the 21 - 
direction, corresponding to the curl of the magnetic field, 
such that Ampere’s law is satisfied 

jz = ~ !)• (2) 

Here, a greater or smaller than unity corresponds to 
a contraction along the x or y axis respectively, and re¬ 
sults in an inwards j x B force on the plasma along the 
same axis. This in turn pushes the field lines inwards, 
which serves to increase the initial imbalance, which in 
turn increases the inwards force and the field collapses. 
Due to the frozen-in condition, this leads to a build up of 
plasma near the X-point and eventually to the formation 
of a diffusion region and a current sheet, where field lines 
reconnect. Since the simulation allows for kinetic effects, 
ions decouple, at a typical length scale of the ion inertial 
length (c/ujpi), allowing further compression of the field 
near the X-point, allowing for fast reconnection. The 
width of the current sheet is given by the electron iner¬ 
tial length, c/ujpe , while the length of the current sheet, 
which is a decisive factor in the reconnection rate, was 
shown to be of the order of the length of the ion diffusion 
region for tearing-mode reconnection f2~il . The dominant 
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term in the generalized Ohm’s law, which results in the 
breaking of the frozen in condition, was shown to be the 
off-diagonal terms of the electron pressure tensor diver¬ 
gence, due to electron meandering motion (see Ref. 0 )- 
This makes the reconnection process fast compared to 
the resistive MHD, which is too inefficient (see Ref. 23|, 
chapter 7.1.1). 

In all our simulation runs, a was set to 1.2 and Bq was 
adjusted such that the Alfven speed at the y-boundary 
was fixed as V a = 0.1c, where c is the speed of light. Fur¬ 
ther, the parameters in the simulation were adjusted to 
assure that the electron plasma frequency at the bound¬ 
ary was equal to the electron cyclotron frequency at the 
boundary, i.e. w pe = w ce . The number density of both 
electrons and ions in the simulation domain, n e and n p , 
was set to 10 16 m~ 2 , while the temperature for both elec¬ 
trons and ions, T e and T p , was set to 6.0 x 10 7 K, matching 
conditions of flaring in the solar corona. The code was 
set to use 200 particles per species per cell, which was 
shown to be sufficient for accurate electromagnetic field 
dynamics in convergence tests. 

The simulation code used here is a 2.5D relativistic and 
fully electromagnetic Particle In Cell (PIC) code, devel¬ 
oped by the EPOCH collaboration, based on the original 
PSC code by Hartmut Ruhl j25|. This code was modi¬ 
fied to allow for closed boundary conditions that conserve 
electromagnetic flux at the boundary. To achieve this, 
zero-gradient boundary conditions are imposed both on 
the electric and magnetic fields in x- and y-directions and 
the tangential component of electric field was forced to 
zero, while the normal component of the magnetic field 
was kept constant. This anchors magnetic field lines on 
the boundary, thus inhibiting loss (or gain) of magnetic 
flux from the simulation domain. This is applicable to 
magnetic fields which are anchored in the solar photo¬ 
sphere by the frozen-in condition and form an X-point 
higher up in the corona, which serves as a simplest model 
of a solar coronal active region. All numerical runs in this 
study use closed boundary conditions (simulation results 
for X-point with open boundary conditions can be found 

in Refs. 17(1 [2~iT|. 


III. OCTUPOLAR STRUCTURE FOR 
DIFFERENT DOMAIN SIZES 

While having the advantage of simulating a self- 
contained, energy conserving system, flux conserving 
boundary conditions naturally have a limiting effect on 
the movement of field lines, especially closer to the edge 
of the domain. In order to diminish this effect and allow 
natural reconnection dynamics to occur, different system 
sizes were used for the square simulation domain, ranging 
from 4 c/cjpi to 32 c/cj p i. The effective grid sizes used 
ranged from 2 L = 2.14 m to 2 L = 17.12 m. For most 
of the runs, simulation grid cells were set to the Debye 
length (Ad), thus leading to grids ranging from 400x400 
to 1600x1600 grid cells. However, for the case of a system 


size of 32c/w p i this was not computationally possible. An 
under-resolved simulation run, using 1600x1600 grid cells 
and thus cell size of 2 Ad, however proved consistent and 
was shown to be energy conserving despite possible nu¬ 
merical heating and is thus included here. The ion mass 
was set to 100 times the electron mass, i.e. rrii = 100 m e , 
to speed up the code. The value of B 0 for the differ¬ 
ent runs was adjusted such that the Alfven speed at 
the y-boundary was fixed as v a = B^ / yj ^orriirii = 0.1c, 
where B^ represents the strength of the magnetic field at 
{pcmaxi 0)- For meaningful comparison between the runs, 
the a;-axes of plots showing time dynamics use Alfven 
time, 

t a = L/v a , (3) 

where L denotes the half-length of the system size. 



FIG. 2: The reconnected flux for simulation runs of different 
domain sizes, given by the difference between magnetic flux 
at X-points and O-points, indicative of the reconnection rate 
(shown for different electron to ion mass ratios in Fig. 2 of 
Ref. [l9}). The solid, dashed, dash-dotted and dotted curves 
show the reconnection rate for domain sizes of 4c/ai p i, 8c/ui p i, 
IQc/ujpi and 32c/w p i respectively. The reconnected flux is 
normalised by Bbc/uipi. For all system sizes a similar amount 
of reconnected flux is reached within 1.5t a . Note that, for 
each case, initial reconnection rate maxima are reached also 
within this period. 

By plotting the reconnected in-plane magnetic flux 
over time for runs with different system sizes (see Fig. 0, 
it was determined that in all runs, a similar amount of 
flux is reconnected within 1.5f a . For all cases an instant 
of maximum reconnection rate is reached within 1.5t a , 
which occurs at progressively later times for greater sys¬ 
tem sizes (note that time normalisation is as stated in Eq. 
©)■ This can be understood as a result of the increas¬ 
ing difference in length between the system size and the 
diffusion region. I.e. while the system size increases, the 
width of the diffusion region remains the same and field 
lines have to travel a greater distance before undergoing 
reconnection. For a system size of 2 L = 4c/w P i the the 
peak reconnection rate is lower than for the other cases, 
for which the peak reconnection rate is approximately 
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FIG. 3: The magnetic out-of-plane field for 
simulation runs with different domain sizes 
at peak reconnection rate, as indicated in 
Fig. m Panels (a) to (d) show runs using 
system sizes of 4 c/cj p i, 8c/uj p i, lGc/ujpi and 
32c/u .’ p i respectively. White lines superim¬ 
posed on plots show the in-plane magnetic 
field at the same time. Different times cor¬ 
respond to time instants when maximum re¬ 
connection rate is reached. Note that Bb is 
the same as B j,. 


equal. In this case, the closeness of the boundary to the 
diffusion region is limiting the reconnection. 

The out-of-plane magnetic field for different system 
sizes at peak reconnection rates is plotted in Fig. [3] The 
width of the out-of-plane magnetic structure does not 
change with the size of the domain. For system sizes of 
8 c/uipi, and greater, it can be seen that the horizontal ex¬ 
tent of the out-of-plane field is approximately contained 
within —4c/w P i to +4 c/ui p i. Since the particle density 
and mass ratios are fixed for the different runs, the width 
of ion diffusion region, ~ c/w P j, is also fixed. Since ions 
and electrons only move independently within the ion dif¬ 
fusion region, which allows for in-plane currents and thus 
out-of-plane magnetic fields to be generated, it is evident 
that the emergence of octupolar structure is an aspect of 
ion diffusion region physics. 

The vertical extent of the out-of-plane magnetic field 
however does increase for greater system sizes. This 
shows that the generation of in-plane currents occurs at 
a vertical distance determined by the strength of the in¬ 
plane magnetic field, which has equal strengths on rela¬ 
tive positions on the domain. The in-plane field lines in 
Fig. [3] show that the value of B x in the outflow region 
is fixed at approximately halfway between the X-point 
and the system boundary, for all simulation runs. Since 
flux can not escape from the boundary, this shows that 
the movement of field lines is halted at the same rela¬ 
tive position on the domain for all system sizes. Thus, 


the motion of electrons must also be halted at this point 
since they recouple to the field shortly after reconnec¬ 
tion. However, ions are still decoupled at this point and 
over-shoot, thus creating in-plane currents and charge 
separation. In response, electrons move along the field- 
lines to restore charge neutrality, resulting in the cur¬ 
rent loops and out-of-plane magnetic field as illustrated 
in Fig. [T] (right). In panel (a) of Fig. [3] we see that for 
a case of a simulation size of 2 L = 4c/uj p i , the octupo¬ 
lar components are not contained within the domain. It 
is thus clear that a simulation size of at least 8 c/co P i is 
required to comprehensively model the dynamics within 
the ion diffusion region, without spurious influence from 
the boundary condition. 


IV. THE OCTUPOLAR FIELD STRUCTURE 
FOR DIFFERENT ION MASSES 

In order to determine that the results relating to oc¬ 
tupolar structure hold true for realistic electron-to-ion 
mass ratios, the relationship between the strength of the 
octupolar and quadrupolar components and the mass ra¬ 
tio was investigated. In these simulation runs, the mass 
of electrons was fixed as their physical mass, m e , and 
ion masses were varied from rrii = 50m e to nii = 400m e . 
In order to simulate comparable effects, runs with dif¬ 
ferent ion masses were adjusted in scale, such that the 
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domain size remained the same in number of ion iner¬ 
tial lengths, c/ujpi = c\Jmi€o/riie 2 . Thus, for greater ion 
masses, the domain size was increased accor ding ly (see 
other findings relating to these setups in Ref. [l9|]). Fur¬ 
ther, as before, the Alfven speed at the boundary was 
fixed (v a = Iyj = 0.1c), which means that B b 
was adjusted accordingly, while Bq remained fixed. For 
meaningful comparison between the runs, time normali¬ 
sations as in Eq. m are applied. 
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FIG. 4: (Top) The strength of the quadrupolar field com¬ 
ponents for runs of different ion masses, and corresponding 
setups as described in the text. Different line styles cor¬ 
respond to different mass ratios as indicated. (Bottom) as 
above, showing the strength of octupolar field components. 

By running the simulations described above, for a 
system size of 4c/w P j, plots in Fig. [I] were obtained. 
Here, the strength of the qudarupolar components 
was determined as the maximum field strength of the 
quadrupolar field at each time step, and similarly for 
the octupolar components, i.e. max(\B z , qua d(x,y)\) and 
m&x(\B ZiOCt (x,y)\) respectively. We gather from Fig. 0 
that increasing the ion mass results in a minor increase in 
the peak strength of the quadrupolar and octupolar field 
components. In all cases, the octupolar components con¬ 
stitute a significant fraction (i.e. 10—15%) of the in-plane 
magnetic held at the boundary and should therefore be 
of significance for observations in laboratory plasma ex¬ 
periments and spacecraft observations. 


V. OCTUPOLAR STRUCTURE FOR 
TEARING-MODE VS X-POINT COLLAPSE 

The primary factor in the emergence of octupolar 
structure in X-point collapse is the role of ion dynamics. 
Due to the large mass of ions compared to electrons, their 


contribution to the out-of-plane Held has generally been 
considered negligible. However, by revisiting tearing¬ 
mode reconnection studies, exhibiting quadrupolar struc¬ 
ture, it was possible for us to show that a similar octupo¬ 
lar structure is also present in these scenarios. Using the 
reconnection setup from kinetic simulation by Pritchett 
i for a tearing-mode reconnection, and carefully inves¬ 
tigating the out-of-plane magnetic Held structure, it was 
shown that octupolar components also emerge here (see 
Fig. 0. 

The initial reconnection magnetic field used in this sim¬ 
ulation is a Harris neutral sheet configuration, given by 

B x = B(,tanh(y/w), (4) 

together with ion and electron density profiles of 

n e ,i = n' 0 sech 2 {y/uj) + n b , (5) 

where B' 0 is determined by the Alfven speed, which is set 
as v a = c/20. The half-thickness of the current sheet, w, 
is set to 0.5 c/uipi, n' 0 is determined by the equilibrium 
condition for the neutral sheet and n b is a constant back¬ 
ground density of 0.2 ?Tq. The mass of ions in the simula¬ 
tion is set as mi = 25m e . By introducing an initial flux 
perturbation, a magnetic island with a transverse size 
comparable to u> is generated, leading to a reconnection 
region comparable to the size of the domain. 

As shown in Fig. 0 there are distinct differences in 
the out-of-plane magnetic field structure for the two re¬ 
connection setups, which could lead to different observa¬ 
tional signatures. In X-point collapse, the quadrupolar 
field emerges in the centre of the domain and two sets 
of octupolar field components emerge in each of the out¬ 
flow and inflow regions. However, in the tearing-mode 
case, a quadrupolar field forms at the centre, but octupo¬ 
lar components emerge in the outflow region only. Using 
line plots, the out-of-plane magnetic field along tracks 
through the octupolar magnetic field structure of the 
X-point collapse and tearing-mode scenarios is shown. 
The selected tracks, in all cases, show two consecutive 
troughs and peaks in the out-of-plane field. In X-point 
collapse, for both the horizontal and vertical track, the 
initial through and the final peak have a smaller field 
strength than the intermediate ones. This is the result of 
the track starting and ending in the outer region of the 
domain, where octupolar components, which are lower 
in field strength, dominate. In the inner region, the out- 
of-plane magnetic field is dominated by the quadrupolar 
components and greater field strengths are observed. For 
the tearing-mode case, for the horizontal track, this trend 
is also observed, but for the vertical track it is reversed, 
i.e. the initial trough and the final peak have a greater 
field strength than the intermediate ones. This can also 
be shown to be the result of the order in which the track 
passes quadrupolar and octupolar field components. In 
the latter case, the track passes through the quadrupolar 
components before and after the octupolar components, 
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FIG. 5: (Left panels) The out-of-plane mag¬ 
netic field structure at peak reconnection for 
X-point collapse (t/t a = 1.1) and tearing 
mode (tui' ci = 20), as in Ref. ||. In both 
cases, elements of opposite polarity can be 
seen next to the inner quadrupolar structure, 
making an overall octupolar field. (Right 
panels) The line-profile of the out-of-plane 
magnetic field at the dotted tracks shown in 
the corresponding contour plots. These pro¬ 
files effectively represent possible observa¬ 
tions by a spacecraft mission, passing though 
reconnection regions of X-point collapse or 
tearing-mode type. 


thus leading to lower field strength in the middle of the 
track. If a spacecraft mission, such as that discussed in 
Ref. (HI , were to pass through a reconnection region and 
observe one of these line profiles, it would be possible to 
distinguish between the possible reconnection setups ac¬ 
cordingly. 

As in Ref. [HI , to determine the causes of the observed 
magnetic poles in terms of the currents in the simulation, 
Ampere’s law was taken in component form, such that 


was possible to individually integrate over the simula¬ 
tion grid, using the three different currents, to obtain 
their individual contributions to B z . E.g. using Eq. © 
one obtains 

r L y 

Bz,ion(Qi Ly) — / y) dy (8) 

Jo 

followed by 0 to get 


dB z = HojxAondy + Hojx,electrondy + ——^-dy (6) 

c- at 

and 

1 dE 

dB z fJ , ojy,iondx f^ojy,electrondx ^ dx. (T) 

By integrating over the domain for these terms, con¬ 
tributions to B z from individual currents, i.e. electron, 
ion and displacement current, can be calculated. The 
individual currents for electrons and ions are provided 
by the simulation at each grid cell and the displacement 
current was obtained by taking a five-point stencil using 
electric field values at the same cell, separated over four 
time steps. As a starting point for the integration a neu¬ 
tral point in B z had to be chosen. For X-point collapse, 
the most logical point was the centre of the domain, since 
here B z (0 , 0) = 0.0. For the tearing-mode case, any point 
point far out in the inflow region was suitable, as these 
regions are shown to be far and disconnected from the 
diffusion region (i.e. field lines from outside the diffu¬ 
sion region here do not enter the diffusion region) and 
thus lacking out-of-plane magnetic structure. Thus, it 


B z ,ion{Lx, Ly) — Rjz,ion(0, Ay) / Jy ,ion (•£; Ay ) dx, 

Jo 

(9) 

where (L x ,L y ) represents an arbitrary point on the B z 
grid. 

Carrying out this integration for all L x and L y on the 
grid for all current contributions for both A'-point col¬ 
lapse and tearing mode, plots shown in Fig. [6] are ob¬ 
tained. For considerations of symmetry, only the lower 
left quarter of the simulation domain is shown (also note 
that, since X-point collapse, by convention, has inflow 
regions along the horizontal axis and tearing-mode along 
the vertical axis, the sign of the quadrupolar and oc¬ 
tupolar field structure is reversed). Top and middle pan¬ 
els show ion and electron current contributions to the 
out-of-plane magnetic field and the bottom panel the 
combined contributions, which lead to an out-of plane 
field as shown in Fig. [5] As expected, based on Hall dy¬ 
namics, the contribution to the quadrupolar components 
is provided entirely by the electron currents in bothre- 
connection scenarios. Also, as established in Ref. [n|, 
the octupolar components are the result of ion currents. 
However, unlike in A'-point collapse, in the tearing mode 
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FIG. 6: (Left three vertical panels) The electron and ion cur¬ 
rent contributions to the out-of-plane magnetic field for an 
X-point collapse scenario with a domain size of 8 c/uj p i at 
peak reconnection rate (t/t a = 1.1), calculated from electron 
and ion currents respectively, based on Ampere’s law. Super¬ 
imposed arrows indicate the strength and direction of currents 
in the domain. Arrows next to plots show current strengths 
corresponding to charged particles of the characteristic den¬ 
sity, moving at the Alfven speed. (Right three vertical panels) 
The same for the tearing-mode setup (tu>' ci = 20) according to 
Ref. 3]. Since ion and electron contributions strongly cancel, 
i.e. were much greater than the resulting field, different scales 
are used as indicated. 


scenario, octupolar components only emerge in the out¬ 
flow region. Further, in the tearing mode case, due to 
a lack of asymmetry in the inflow region, there are two 
great contributions from both the ion and electron cur¬ 
rents, which cancel to give a neutral field at the edge of 
the domain. 

Ref. [23|], chapter 5, gives a detailed analysis of the 
differences in reconnection setups with almost uniform, 
i.e. straight inflowing field lines, and non-uniform, i.e. 
curved field lines, reconnection scenarios. From the re¬ 
sults shown in this study, we find that differences in these 
scenarios can be extended to out-of-plane magnetic field 
structure. From Fig. [6]it can be seen that, for the tearing¬ 
mode scenario, there is a uniform current inflow of both 
ions and electrons in the inflow region, which is the rea¬ 
son no out-of-plane magnetic structure is generated in 
this region. This uniformity in particle inflow is a di¬ 
rect consequence of the lack of curvature of the inflow¬ 
ing field lines. In contrast to this, in X-point collapse, 



FIG. 7: The setup and observed resulting current and out-of- 
plane magnetic field generation for a tearing-mode reconnec¬ 
tion scenario. As in Fig.[l]labelled tracks illustrate the motion 
of electrons and ions, generating the out-of-plane magnetic 
field, and black arrows on field lines indicate the inflow and 
outflow directions. 


the the curvature in the field-lines means that paths of 
ions and electrons greatly diverge when ions decouple and 
thus currents and out-of-plane fields, as shown Fig. [[) 
are generated. Further, from Fig. [Gj we can see that 
electrons in the tearing-mode scenario exhibit the same 
Hall-dynamics and flow along the field-lines towards and 
away from the X-point when entering the diffusion re¬ 
gion. However, the field-lines in the outflow region are 
no longer straight and electrons move in a curved path, 
away from the X-line, while decoupled ions move in a 
straight horizontal path away from the X-point. It is this 
divergence in ion and electron flows which leads to the 
current loops that generate the octupolar components in 
the out-of-plane magnetic field structure in the tearing¬ 
mode case. The difference of this generation mechanism 
to that shown for A'-point collapse (see Fig. [T]) can be 
seen in Fig. [7] 


VI. CONCLUSIONS 

Results relating to the generation of out-of-plane oc¬ 
tupolar structure in collisionless reconnection in an X- 
point collapse [16] were extended by the investigation of 
the effect of changing domain size and ion masses. It 
was established that, when fixing the Alfven speed at 
the boundary and increasing the size of the domain, the 
horizontal extent of the octupolar region remains approx¬ 
imately the same (see Fig. O. This was found to be con¬ 
sistent with previous findings since octupolar structure 
was shown to be the result of in-plane ion currents. As 
ion currents can only contribute to the out-of-plane mag¬ 
netic field within the ion diffusion region and, since c/co p i 
remained fixed for different system sizes, the width of the 
octupolar region remained fixed as well. The length of 
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the octupolar structure however increased for greater sys¬ 
tem sizes. From the in-plane field in Fig. [3]it was possible 
to see that the outward motion of field lines in the out¬ 
flow region diminishes at the same relative position on 
the domain for all system sizes. As ions remain decou¬ 
pled at this point, they over-shoot, causing in-plane cur¬ 
rents. Further, electrons compensate for this by moving 
along the field lines, thus resulting in the current loops 
and out-of plane magnetic field as shown in Fig. |Tj 

By varying the electron to ion mass ratio and adjusting 
the simulation domain size to keep the same number of 
ion inertial lengths, it was shown that octupolar field 
strength was not significantly affected by the electron 
to ion mass ratio variation and consistently remained a 
significant fraction of the in-plane field, i.e. between 10% 
and 15% (see Fig. |lj). Hence, octupolar structure in an 
X-point collapse would also be generated in plasmas with 
realistic mass-ratios. 

Further, the discovery of octupolar structure in a 
tearing-mode scenario is presented and the differences in 
the generation process are analysed. It is found that, due 
to the uniform nature (i. e. straightness of field lines) of 
tearing-mode collapse, significant octupolar structure is 
only generated in the outflow region (see Fig. ED. Here, 
ions flow out along the centre of the diffusion region, 
while electrons follow a curved path along the field lines 
(as shown in Fig. ED and thus in-plane currents and out- 
of-plane fields are formed. Further, by analysing the line 
profiles of the out-of-plane magnetic field along potential 
tracks through the reconnection region, distinctly differ¬ 
ent profiles were obtained which could have relevance to 
the identification of magnetic structures found by space¬ 


craft missions (see Fig. O. 

It is shown in Ref. [26j, Fig. 7, that plasma flow along 
the current sheet can induce streaming sausage and kink 
modes that generate out-of-plane magnetic fields resem¬ 
bling the one discussed here. Further, in Ref. f27j a 
framework is presented where the Hall fields of recon¬ 
nection are formulated as a process of Alfven eigenmode 
generation and dissipation. It was shown that for the 
n = 1 mode that this frame work provides analytical pre¬ 
dictions for the Hall fields, which are in good agreement 
with observed magnetic fields. In higher wave modes 
(n>l) this framework returns a higher order structure 
in the Hall fields which could potentially be the cause of 
octupolar structure in the out-of-plane magnetic field ob¬ 
served in this study. The application of this framework in 
the X-point collapse case is currently under investigation 
and will be reported elsewhere. 
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